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and search for a Dark Matter signal 
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Abstract 



Observations of the Sagittarius dwarf spheroidal (Sgr dSph) galaxy were carried out 
with the H.E.S.S. array of four imaging air Cherenkov telescopes in June 2006. A 
total of 11 hours of high quality data are available after data selection. There is 
no evidence for a very high energy 7-ray signal above the energy threshold at the 
target position. A 95% C.L. flux limit of 3.6 x 10 _12 cm _2 s _1 above 250 GeV has 
been derived. Constraints on the velocity-weighted cross section (uv) are calculated 
in the framework of Dark Matter particle annihilation using realistic models for the 
Dark Matter halo profile of Sagittarius dwarf galaxy. Two different models have 
been investigated encompassing a large class of halo types. A 95% C.L. exclusion 
limit on (av) of the order of 2 x 10~ 25 cm 3 s _1 is obtained for a core profile in the 
100 GeV - 1 TeV neutralino mass range. 

Key words: Gamma-rays : observations - Dwarf Spheroidal galaxy, Dark Matter 
PACS : 98.70.Rz, 98.56.Wm, 95.35.+d 
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1 Introduction 



Astrophysical and cosmological observations provide a substantial body of ev- 
idences for the existence of Cold Dark Matter (CDM) although its nature 
remains still unknown. It is commonly assumed that CDM is composed of yet 
undiscovered non-baryonic particles for which plausible candidates are Weakly 
Interacting Massive Particles (WIMPs). In most theories, candidates for CDM 
are predicted in theories beyond the Standard Model of particle physics [1]. 
The indirect detection of Dark Matter (DM) annihilation may bring new in- 
sights to probe the astrophysical nature of Dark Matter. To a large extent 
complementary to direct searches, the indirect detection enables to search for 
DM outside the solar system. Indeed, it may give detailed morphology features 
that may constrain the DM halo profile. The annihilation of WIMPs into 7- 
rays may lead to detectable very high energy (VHE, E > 100 GeV) 7-ray 
fluxes above background via continuum emission from hadronization of gauge 
bosons and heavy quarks, or 7-ray lines through loop-induced processes. The 
H.E.S.S. (High Energy Stereoscopic System) array of Imaging Atmospheric 
Cherenkov Telescopes (IACTs), designed for high sensitivity measurements in 
the 100 GeV - 10 TeV energy regime, is a suitable instrument to detect VHE 
7-rays and investigate their possible origin. 

Various astrophysical systems ranging from local objects in the galactic halo to 
galaxy cluster scales have been considered as targets for DM annihilation 7-ray 
studies. Interacting with baryonic matter only through gravity, the WIMPs 
are expected to concentrate at the centre of high density region. A probable 
candidate is the elliptical galaxy M87 at the centre of the Virgo cluster [2]. 
However, the temporal variability of the H.E.S.S. signal [3] excludes the bulk 
of the signal in the TeV range to be of a dark matter origin. Prospects of indi- 
rect detection from the Andromeda galaxy M31 [4] and the Large Magellanic 
Cloud [5] have been also investigated. H.E.S.S. observations in the direction 
of the Galactic Center (GC) have revealed a source of VHE 7-ray emission 
(HESS J1745-290). The measured spectrum is difficult to reconcile with a 
DM interpretation and does not match the expected annihilation spectrum. 
The very high energy cut-off above 7 TeV requires an uncomfortably massive 
DM particle [6]. Moreover, standard astrophysical objects found in the region 
such as the supermassive black hole Sgr A* or the recently discovered ple- 
rion G359.95-0.04, can easily account for the observed signal. The shell-type 
supernova remnant Sgr A East is unlikely to be associated to the signal [7]. 
Besides HESS J1745-290, the deep observations carried out with H.E.S.S. in 
the GC region have highlighted the existence of a diffuse component along the 
galactic ridge [8]. This extended TeV emission may be explained by cosmic 
ray interactions inside the central molecular zone leading to diffuse astrophys- 
ical backgrounds that might hide the exotic signal [9]. Future indirect DM 
searches in this region will have to overcome this challenging background. In 
contrast, dwarf galaxies may be a less complex environment because of the 
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reduced amount of gas in such small systems [10]. 

Dwarf Spheroidal galaxies (dSph) such as Sagittarius dwarf or Canis Major, 
discovered recently in the Local Group, are among the most extreme DM- 
dominated environments. Indeed, measurements of roughly constant radial 
velocity dispersion of stars imply large mass to luminosity ratios [11]. Nearby 
dwarfs are ideal astrophysical probes of the nature of DM as they usually con- 
sist of a stellar population with no hot or warm gas, no cosmic ray population 
and little dust. The Sagittarius dwarf Spheroidal galaxy (Sgr dSph) is the 
next-to-last Galactic satellite galaxy discovered [12]. The core of Sgr dSph is 
located at 1 = 5.6° and b = -14° in galactic coordinates at a distance of about 
24 kpc from the Sun [13]. Sgr dSph has made at least ten Milky Way cross- 
ings it should thus contain a substantial amount of DM to avoid to have been 
entirely disrupted. Latest velocity dispersion measurements on M giant stars 
with 2MASS yields a light to mass ratio of about 25 [14]. The Sgr dSph core is 
positioned behind the bulge of Milky Way but outside the Galactic plane, thus 
reduced foreground 7-ray contaminations are expected. The luminous density 
profile of Sgr dSph has two components [15]. The compact component, namely 
the core, is characterized by a size of about 3 pc FWHM, which corresponds 
to a point-like region for H.E.S.S. This is the DM annihilation region from 
which 7-ray signal may be expected. A diffuse component is well fitted by a 
King model with a characteristic size of 1.6 kpc. 

In this paper, we present the observation of the Sagittarius dwarf galaxy by 
the H.E.S.S. array of Imaging Atmospheric Cherenkov Telescopes based on a 
dataset collected in June 2006. A careful modeling of the Dark Matter halo us- 
ing the latest measurements on the structural parameters of Sagittarius dwarf 
is presented to derive constraints on the WIMP velocity-weighted annihilation 
rate. 



2 Search of VHE 7-rays from observations of Sagittarius Dwarf by 
H.E.S.S. 

2.1 The H.E.S.S. instrument 

H.E.S.S. (High Energy Stereoscopic System) is an array of four imaging at- 
mospheric Cherenkov telescopes located in the Khomas Highlands of Namibia 
at an altitude of 1800 m above sea level. Each telescope consists of an optical 
reflector of about 107 m 2 effective area composed of 382 round mirrors ar- 
ranged on a Davis-Cotton mount [16]. The interaction of the primary 7-ray in 
the Earth's upper atmosphere initiates an electromagnetic shower. The reflec- 
tor collects the Cherenkov light emitted by the charged particles composing 
this shower, and focuses it onto a camera comprising 960 fast photomultipli- 
ers (PMTs) of individual field of view of 0.16° diameter [17]. Each tube is 
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equipped with Winston cones to limit the field of view of each PMT and min- 
imize the background light. The total field of view of the H.E.S.S. instrument 
is 5° in diameter. The stereoscopy technique used in the imaging atmospheric 
Cherenkov telescopes allows for accurate reconstruction of the direction and 
energy of the primary 7-rays as well as an efficient rejection of the background 
induced by cosmic ray interactions [18]. The energy threshold of H.E.S.S. at 
zenith before selection cuts moved from 100 GeV at the commissioning of the 
experiment in 2003 to 160 GeV due to the degradation of the optical perfor- 
mance in 2006. The point source sensitivity is better than 2 x 10 _13 cm _2 s _1 
above 1 TeV for a 5a detection in 25 hours [19]. 



2.2 Dataset 



The observations of the Sgr dSph were taken in June 2006 with zenith angles 
ranging from 7° to 43° around an average value of 19°. Data were taken in 
28-minute observation runs in the wobble mode method with pointing direc- 
tions offset by an angular distance of typically ±0.7° from the nominal target 
position. The dataset suitable for analysis was selected using the standard run 
selection procedure [19], which in particular removes the data taken under bad 
atmospheric conditions. A total of 25 runs out of 26 are selected for the anal- 
ysis. After calibration of the raw shower images from PMT signals [20], two 
independent reconstruction techniques were combined to select 7-ray events 
and reconstruct their direction and energy. The first one uses the Hillas mo- 
ment method [21]. The second analysis referred hereafter as "Model Analysis", 
is based on the pixel-per-pixel comparison of the shower image with a template 
generated by a semi-analytical shower development model. The event recon- 
struction relies on a maximum likelihood method which uses available pixels 
in the camera, without requirement for an image cleaning [22,23]. The recon- 
structed shower parameter (energy, impact, direction and primary interaction 
point) are obtained as a product of the minimization procedure. The sepa- 
ration between 7 candidates and hadrons is done using a combination of the 
Model goodness-of-fit parameter [23] and the Hillas mean scaled width and 
length parameters, which results in an improved background rejection [19]. 
Standard cuts on the width and the length of Hillas ellipses combined with 
the goodness-of-fit are used to suppress the hadronic background [21]. An ad- 
ditional cut on the primary interaction depth is used to improve background 
rejection. Both methods yield a typical energy resolution of 15% above energy 
threshold. In the Model Analysis, the angular resolution at the 68% contain- 
ment radius is found to be better than 0.06° per 7-ray. 



7 



2.3 Data analysis 



The on-source signal is denned by integrating all the events with angular 
position 9 in a circle around the target position with a radius of 9 cut . The target 
position is chosen according to the photometric measurements of the Sgr dSph 
luminous cusp showing that the position of the centre corresponds to the centre 
of the globular cluster M54 [24] . The target position is thus found to be (RA = 
18 h 55 m 59.9s, Dec = -30d28'59.9") in equatorial coordinates (J2000.0) or (1 = 
5°41'12.9",b = -14°16'29.8") in Galactic coordinates. The signal coming from 
Sgr dSph is expected to come from a region of 1.5 pc, about 30", much smaller 
the H.E.S.S. point spread function (PSF). A 9 cut value of 0.14° suitable for a 
point-like source was therefore used in the analysis. In case of a Navarro-Frenk- 
White (NFW) density profile [25] for which p follows r _1 or a cored profile [10] 
folded with the point spread function (PSF) of H.E.S.S., the integration region 
allows to retrieve a significant fraction of the expected signal. See Table 1. A 
cut on the image size of 60 photoelectrons is used to obtain a good sensitivity 
for weak sources. In order to reduce systematic effects which affect images 
close to the edges of the camera, only events reconstructed within a maximum 
distance of 2.5° from the camera centre are used for this analysis. The excess 
sky map is obtained by the subtraction of a background model on the 7- 
ray candidate sky distribution. The background level is estimated using the 
ring-background method [26] where the background rate is calculated from 
the integration of 7-like events falling in an annulus around the centre of the 
camera with identical observation conditions and acceptances than that used 
for the on-source region, which allows an estimate of the background on every 
sky position. 

The excess sky map of the 7-ray candidates is presented in Fig. 1 in the 
RA/Dec J2000 coordinates centered on the Sgr target position. For each bin, 
7-ray like events are summed within a radius of 0.14°. No 7-ray excess is found 
at the target position. Further, no significant excess is observed anywhere else 
in the sky map. The 9 2 distribution of the observed 7-ray events relative to the 
target position, is presented in Fig. 2 as well as the background distribution. 
In order to check the robustness of the results presented here, the dataset 
has been analyzed using different analysis methods as the so-called Hillas 
and model3D [27] methods. No significant 7-ray excess is detected in the 
corresponding sky maps. 

Since Fig. 2 shows no 7-ray excess, we derived the 95% confidence level upper 
limit on the observed number of 7-rays : N^ 5%C,L \ The limit is computed 
knowing the numbers of events above the energy of 250 GeV in the signal 
region Nqn = 437, in the background region Nqff = 4270, and the ratio of 
the off-source time to the on-source time a = 10.1. We use the Feldman & 
Cousins method [28] and obtain : 

N 95%C.L. = 5 g (1) 
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Given the acceptance of the detector for the observations of the dSph Sgr, an 
95% confidence level upper limit on the 7-ray flux is also derived : 

$ 7 (E 7 > 250 GeV) < 3.6 x 1(T 12 cm^s" 1 (95%C.L.) 



3 Predictions of 7-rays from Dark Matter annihilation 



3. 1 Theoretical framework 



The annihilation of DM particles can generate 7-ray fluxes through different 
processes depending on the particle physics scenarios. Generally, WIMP an- 
nihilations will produce a continuum of 7-rays with energies up to the WIMP 
mass issued from the hadronization and decay of the cascading annihilation 
products, predominantly from 7r°'s generated in the quark jets. In the R-parity 
conserving supersymmetric extensions of the Standard Model, the lightest su- 
persymmetric particle (LSP) is a stable particle and is a good CDM candidate. 
The LSP is, in various SUSY breaking models, the lightest neutralino % Being 
electrically neutral and colorless, it is among the best motivated candidates to 
account for CDM [29]. In Minimal Supersymmetric Standard Model (MSSM) 
scenarios [29], the annihilation of neutralinos can, on top of the hadroniza- 
tion continuum, produce also monoenergetic spectral lines of 7-rays result- 
ing from loop-induced annihilation processes such as XX ~^ 77; XX ~* 1%, 
XX lh, even though these are very challenging to detect experimentally 
due to the high suppression of such final states. Beyond the Standard Model, 
plausible candidates are provided by the universal extra dimension (UED) 
theories. In Kaluza-Klein (KK) scenarios with KK-parity conservation, the 
lightest Kaluza-Klein particle (LKP) is stable [30], the best-motivated being 
the first KK mode of the hypercharge gauge boson, In this case, pairs 
annihilate preferentially into charged lepton pairs which radiatively produce 
7 with harder spectra. Cascading decays of qq final states lead to secondary 
7-rays [31]. 

The 7-ray flux from annihilations of DM particles of mass uidm accumulating 
in a spherical DM halo can be expressed in the form : 

«BS) = I H^xJ M 4!l (2) 
dE 1 An m 2 DM dE 1 — J. ' 

y v ' Astrophysics 

Particle Physics 

as a product of a particle physics component with an astrophysics component. 
The particle physics part contains (o~v), the velocity- weighted annihilation 
cross section, and dN^/dE^, the differential 7-ray spectrum summed over the 
whole final states with their corresponding branching ratios. The astrophysical 
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part corresponds to the line-of-sight-integrated squared density of the DM 
distribution J, averaged over the instrument solid angle integration region for 
H.E.S.S. (An = 2 x KT 5 sr) : 

J = f p 2 ( r [s])ds J (AQ) = / PSF*Jdtt (3) 

Ji.o.s Ail Jaq. 

where PSF is the point spread function of H.E.S.S. 



3.2 Modeling the Sagittarius dwarf Dark Matter halo 



The mass distribution of the DM halo of Sgr dwarf has been described by 
plausible models taking into account the best available measurements of the 
Sgr dwarf galactic structure parameters. We have used two widely different 
models. The first has a NFW cusped profile [25] with the mass density given 
by: 

j\ 

pNFw(r) = -7— ^ (4) 

r(r + r s y 

with A the normalization factor and r s the scale radius taken from [10]. Using 
Eq. 3, the value of J obtained with this model is reported in Table 1. 
We have also studied a core-type halo model as in [10] characterized by the 
mass density : 

_ v\ 3r* + r 2 

PcoreKT) - ^ Q ^ + ^ W 

where r c is the core radius and v a a velocity scale. However, we have tried to 
update the v a and r c values which were used in [10]. By inserting in the Jeans 
equation the luminosity profile of the Sgr dwarf core of the form : 

= j-r—ty (6) 

we estimated from the data of reference [24] a = 2.69 ± 0.10 and r c = 1.5 pc. 
Note that the value of r c is only an upper limit. The value of the central ve- 
locity dispersion of Sgr Dwarf is a = 8.2 ± 0.3 kms -1 [32]. We have assumed 
that the velocity dispersion is independent of position. The value of v a is then 
given by v a = \fa a = 13.4 kms" 1 . The cored model gives a very large value of 
J, which is reported in Table 1. The third column of Table 1 gives the amount 
of signal expected in the solid angle integration region AQ = 2 x 10~ 5 sr. 
The value of J in the cored model depends on various parameters such as the 
radial dispersion velocity, the baryon fraction in the core, the core radius and 
the velocity tensor anisotropy. The central value of the radial dispersion ve- 
locity was taken as 8.2 kms -1 . Previous measurements report 11.4 kms -1 [33] 
which would lead to an increase of a factor 4 in J. Deviations from the asymp- 
totically flat rotation curve have been studied. The effect of the a parameter 
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from [10] leads to a 50% increase or decrease in J, according to the sign of a. 
The anisotropy in the velocity dispersion may lead to a reduction of a factor 
2 in the dark matter density in the central region with respect to the case 
without anisotropy [34]. The baryon mass fraction in the very central region 
can not be neglected. However, the effect on J turns out to be at most 50% 
due to the relatively fast increase of the mass-to-luminosity ratio with radius. 
Finally, in our model, the ratio of the luminous core radius to the dark matter 
core radius is p = 1. In general, as emphasized in [35], p could be lower than 
1. J is strongly dependent on the value of p. Extending the dark halo to 200 
pc would lower J by 2 orders of magnitude. Note also that our value of the 
luminous radius is only an upper limit. Decreasing the luminous radius leads 
to an increase of J. 



3. 3 Sensitivity 



For a given DM halo, the relevant quantities for DM particle annihilation 
searches are the DM particle mass m DM and the velocity- weighted cross sec- 
tion (crv). With the limit on the number of 7, N 7 , derived in section 2.4, we 
can compute the limit on (crv) from H.E.S.S. results with the Sgr dwarf DM 
halo profile modeled in section 3.2. N 7 may be computed using the formula : 



T oi . f" M A^^^ldE., (7) 



■"1 



where A e g corresponds to the effective area of the instrument obtained from 
Monte Carlo simulations as a function of the zenith angle, the offset of the 
source from the pointing direction, the energy of the event and the selection 
cuts. Using the expression of the differential flux given in Eq. 2 with Eq. 7 
yields a 95% C.L. exclusion limit on the velocity-weighted cross section versus 
the DM particle mass for a given halo profile as defined by : 

4tt m 2 N 95%C.L. 
1 \ 95% C.L. _ 2_ _ m DM (o\ 

Jo Acs{E ^dE; d ^ 

where dN^/dE^ is computed with a parametrization of the differential contin- 
uum photon spectrum from [36] for a higgsino-type neutralino. Fig. 3 shows 
the limits in the case of a cored (green dashed line) and cusped NFW (red 
dotted line) profile using the value of J computed in section 3.2. Predictions 
for phenomenological MSSM (pMSSM) models are displayed (grey points) as 
well as those satisfying in addition the WMAP constraints on the CDM relic 
density ficDMh 2 denoted as blue points. The values allowed by WMAP are 
taken to lie in the range 0.09 < f2 CDM h 2 < 0.11. The SUSY models are cal- 
culated with DarkSUSY4.1 [37] in pMSSM framework and characterized by 
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a basic set of independent parameters : the higgsino mass parameter //, the 
gaugino mass parameter M 2 , the CP-odd Higgs mass M^, the common scalar 
mass m , the trilinear couplings A tjb and the Higgs vacuum expectation value 
ratio tan/5. The set of parameters for a given model is randomly chosen in a 
parameter region encompassing a large class of pMSSM models, as described 
in Tab. 2. 

In the case of a cusped NFW profile, the H.E.S.S. observations do not establish 
severe constraints on the velocity-weighted cross section. For a cored profile, 
due to a higher central density, stronger constraints are derived and some 
pMSSM models can be excluded in the upper part of the pMSSM scanned 
region. 

In the case of KK dark matter, the differential 7 spectrum is parametrized 
using Pythia [38] simulations and branching ratios from [30]. Predictions for 
the velocity-weighted cross section of B^ dark matter particle are performed 
using the formula given in [39]. In this case, the expression for (av) depends 
analytically on the mass square. Fig. 4 shows the sensitivity of H.E.S.S. 
in the case of Kaluza-Klein models where the hypercharge boson BW is the 
LKP, for a cored (green solid line) and a cusped NFW (red solid line) profile 
respectively using the value of J computed in section 3.2. With a NFW pro- 
file, no Kaluza-Klein models can be tested. In the case of a cored model, some 
models providing a LKP relic density compatible with WMAP constraints can 
be excluded. From the sensitivity of H.E.S.S., we derive a lower limit on the 
B« mass of 500 GeV. 



4 Discussion 

The Sgr dSph galaxy is among the best target to search for DM signal. Sgr 
dSph region is devoid of astrophysical background unlike the Galactic Center 
region. Indeed, the presence of several 7-ray emitters and diffuse emission 
makes difficult to disentangle the emission from the very center from that of 
other objects. The absence of gas in Sgr offers a cleaner environment to search 
for 7-ray emission. 

The DM profile uncertainties for Sgr dSph are about one order of magnitude. 
In contrast, the Milky Way DM profile suffers from large uncertainties, up 
to five orders of magnitude, due to the difficulty of measuring the structural 
parameters of our Galaxy in the central 10~ 2 pc. 

A modest observation time allows us for the first time to test some pMSSM 
models for the 7-ray annihilation of neutralinos in the direction of Sgr dSph. 
Some KK models can already be excluded. 

Sgr dSph is a target for deeper observations by H.E.S.S. With 100 hours, 
H.E.S.S. will be able to test pMSSM models assuming a NFW profile. With 
such an observation time, H.E.S.S. could exclude all the cosmologically allowed 
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KK models in case of a core profile. 



5 Conclusion 

The observations of the Sagittarius dwarf spheroidal galaxy with H.E.S.S. 
performed in June 2006 reveal no significant 7-ray excess at the nominal target 
position. The Sagittarius dwarf dark matter halo has been modeled using latest 
measurements of its structure parameters. Constraints have been derived on 
the velocity- weighted cross section of the dark matter particle in the framework 
of supersymmetric and Kaluza-Klein models. 
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Fig. 1. Sky map of the 7-ray candidates with an oversampling radius of 0.14°. No 
excess is observed at the target position (RA = 18h54m40s, Dec = — 30d27m05s) 
in equatorial coordinates (J2000) marked with a black triangle. Other spots in the 
field of view are not significant. 
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Fig. 2. 6 2 radial distribution of the ON and OFF events for 7-ray like events from 
the target position (RA = 18h54m40s, Dec = — 30d27m05s) (black dots). Estimated 
background calculated as explained in the text is shown (black triangles). No excess 
is seen at small 9 2 value. 
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Fig. 3. Upper limits at 95% C.L. on (av) versus the DM particle mass in the case of 
a cusped NFW (red dotted line) and a cored (green dashed line) DM halo profiles 
respectively. The pMSSM parameter space was explored with DarkSUSY 4.1 [37], 
each point on the plot corresponding to a specific model (grey point). Amongst 
these models, those satisfying in addition the WMAP constraints on the CDM 
relic density are overlaid as blue square (see text for details). The limits in case of 
neutralino dark matter from pMSSM are derived using the parametrisation from 
reference [36] for a higgsino type neutralino annihilation 7 profiles. 
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Fig. 4. Upper limits at 95% C.L. on (av) versus the DM particle mass in the 
Kaluza-Klein scenarios for a cusped NFW (red dotted line) and a cored (green 
dashed line) DM halo profiles respectively. The blue line corresponds to Kaluza- 
-Klein models [30]. Overlaid (yellow line) are the KK models satisfying WMAP 
constraints on the CDM relic density (see text for details). 
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Table 1 

Structural parameters for a cusped NFW (r s , A) and a cored (r c , v a ) DM halo model, 
respectively. The values of the solid-angle-averaged l.o.s integrated squared DM 
distribution are reported in both cases for the solid angle integration region AO = 
2 x 10~ 5 sr. 



Halo type 


Parameters 


J 

(10 24 GeV 2 cm- 5 ) 


Fraction of signal 
in AO = 2 x 10" 5 sr 


Cusped NFW halo 


r s = 0.2 kpc 
A = 3.3 x 1O 7 M 


2.2 


93.6% 


Cored halo 


r c = 1.5 pc 
v a = 13.4 km s _1 


75.0 


99.9% 



Table 2 

Region of the pMSSM parameter space randomly scanned to generate the models. 
A set of free parameters in the considered ranges corresponds to a pMSSM model. 

pMSSM parameter space region 



100 GeV < n < 30 TeV 
100 GeV < M 2 < 50 TeV 
50 GeV < M A < 10 TeV 
100 GeV < m < 1 TeV 
-3 TeV < A t < 3 TeV 
-3 TeV < A b < 3 TeV 
1.2 < tan/3 < 60 
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List of Figures 

1 Sky map of the 7-ray candidates with an oversampling 
radius of 0.14°. No excess is observed at the target position 
(RA = 18h54m40s, Dec = -30d27m05s) in equatorial 
coordinates (J2000) marked with a black triangle. Other spots 
in the field of view are not significant. 

2 9 2 radial distribution of the ON and OFF events for 7-ray 
like events from the target position (RA = 18h54m40s, Dec = 
— 30d27m05s) (black dots). Estimated background calculated 
as explained in the text is shown (black triangles). No excess 
is seen at small 9 2 value. 

3 Upper limits at 95% C.L. on (<jv) versus the DM particle 
mass in the case of a cusped NFW (red dotted line) and 
a cored (green dashed line) DM halo profiles respectively. 
The pMSSM parameter space was explored with DarkSUSY 
4.1 [37], each point on the plot corresponding to a specific 
model (grey point). Amongst these models, those satisfying 
in addition the WMAP constraints on the CDM relic density 
are overlaid as blue square (see text for details). The limits in 
case of neutralino dark matter from pMSSM are derived using 
the parametrisation from reference [36] for a higgsino type 
neutralino annihilation 7 profiles. 

4 Upper limits at 95% C.L. on (<jv) versus the DM particle 
mass in the Kaluza-Klein scenarios for a cusped NFW 
(red dotted line) and a cored (green dashed line) DM halo 
profiles respectively. The blue line corresponds to Kaluza-Klein 
models [30]. Overlaid (yellow line) are the KK models 
satisfying WMAP constraints on the CDM relic density (see 
text for details). 
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1 Structural parameters for a cusped NFW (r s ,A) and a cored 
(r c , v a ) DM halo model, respectively. The values of the 
solid-angle-averaged l.o.s integrated squared DM distribution 
are reported in both cases for the solid angle integration region 
An = 2 x 10~ 5 sr. 
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Region of the pMSSM parameter space randomly scanned to 
generate the models. A set of free parameters in the considered 
ranges corresponds to a pMSSM model. 
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